PHYSICS OF THE RHYTHMIC APPLAUSE 
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We report on a series of measurements aimed to characterize the development and the dynamics 
of the rhythmic applause in concert halls. Our results demonstrate that while this process shares 
many characteristics of other systems that are known to synchronize, it also has features that are 
unexpected and unaccounted for in many other systems. In particular, we find that the mechanism 
lying at the heart of the synchronization process is the period doubling of the clapping rhythm. The 
characteristic interplay between synchronized and unsynchronized regimes during the applause is 
the result of a frustration in the systems. All results are understandable in the framework of the 
Kuramoto model. 
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I. INTRODUCTION 

Everybody experienced after a good performance the 
audience showing it's appreciation with a thunderous 
synchronized clapping. This is suddenly arising from the 
initially incoherent clapping of individuals and might dis- 
appear and reappear again several times during the ap- 
plause. This synchronization process in the concert hall 
offers a wonderful example of social self-organization, and 
it is believed to be the human scale example of the syn- 
chronization processes observed in numerous systems in 
nature. As detailed in [[l] spontaneous synchronization 
occurs in many biological and sociological systems: fire- 
flies in Southeast Asia synchronize their flashes; crick- 
ets synchronize their chirping; neural cells of the brain 
synchronize voltage fluctuations; pacemaker cells in the 
heart synchronize their fire; womens living together for 
long times find their menstrual cycle synchronized. 

The accepted way of modeling the synchronization pro- 
cess is by considering coupled non-linear oscillators. Un- 
derstanding the synchronization process of coupled oscil- 
lators is an old problem in physics, mathematics and the- 
oretical biology 0. It is believed that the problem goes 
back to Huygens, who realized first the synchronization 
of pendulum clocks hanging together on a wall. While 
the synchronization of identical coupled oscillators inter- 
acting with phase- minimizing interactions is obvious, the 



problem is non-trivial if one considers a population of 
non-identical oscillators. Depending on the strength and 
type of the interaction, and the dispersion of the oscil- 
lators frequencies, synchronization might, or might not 
appear. For biological populations it is crucial to con- 
sider non-identical internal frequencies, since individuals 
are not rigorously equivalent. In this way in order to 
understand synchronization in biological or sociological 
systems we have to consider the more complex problem 
of coupled non- identical oscillators. 

For a theoretical study of coupled oscillator models 
many interaction types were studied. Considering math- 
ematically coupled maps || the problem became inter- 
esting mostly from the view-point of dynamical systems 
and chaos. Since most of the coupling processes in na- 
ture are through pulse- like interactions (firing of neurons, 
flashes of fireflies, clapping etc.) a very realistic pulse- 
like coupling is considered in the integrate and fire type 
models [||| . In these models the firing of an oscillator re- 
sults in a phase-jump of all other oscillators. Because the 
phase-evolution of each oscillator is considered non-linear 
in time, under very general conditions the jumping pro- 
cess leads to synchronization. Interactions continuously 
present in time and thus memory effects are considered in 
the models introduced by Winfree || and Kuramoto . 
In these later models we have phase-difference minimiz- 
ing interactions between globally coupled rotators. For 
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a suitable chosen interaction form, the problem becomes 
exactly solvable and leads to very interesting results. 

Although at a first choice the integrate and fire type 
models look the most promising to understand the syn- 
chronization and dynamics of the rhythmic applause, we 
will argue in favor of the phase-coupled Kuramoto model. 
According to our interpretation of the phenomenon in 
the appearance of rhythmic applause memory effects are 
crucial. By considering the pulse-like interaction of inte- 
grate and fire models this important ingredient is totally 
neglected. In contrast, the continuous phase-coupling in 
the Kuramoto model offers a first approach to deal with 
the relevant memory effects. 

The purpose of our present paper is to reveal some 
interesting and new peculiarities of the rhythmic ap- 
plause, and to understand the observed phenomenon in 
the framework of the well-known Kuramoto model. Of 
course, in the context of the rhythmic applause the Ku- 
ramoto model should be viewed as an abstract theoretical 
representation of the situation, containing the minimum 
number of assumptions already allowing the kind of syn- 
chronization exhibited by the spectators. We mention 
that the main results of our study was briefly discussed 
in a very short recent communication j8|. Here we will 
consider a much more detailed and argumented presen- 
tation. 



II. THE KURAMOTO MODEL 

In order to construct a simple physical model for syn- 
chronization one can take a population of non-identical 
rotators globally coupled through phase-difference min- 
imizing interactions. Synchronization in this system is 
not obvious at all and triggered the interesting studies 
of Winfree [|| . He found that mutual synchronization is 
possible or not depending on the relation between the 
oscillators frequency distribution width (dispersion) and 
the strength of coupling. Later, Kuramoto proposed, and 
together with Nishikava exactly solved Q an elegant re- 
formulation of Winfree's model. 

In the Kuramoto model we have N rotators each of 
them described by it's 4>k phase. The rotators have a 
g(u>) distribution of their lok natural frequencies. Ev- 
ery rotator interacts with all the other ones via phase- 
difference minimizing terms: 
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The N coupled differential equations describing the over- 
damped oscillators dynamics are: 
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Mathematically the synchronization level will be charac- 
terized by an order parameter, q, defined at any time- 
moment as: 



1 



N 

N 

3 = 1 



(3) 



The maximal possible value q — 1 corresponds to total 
synchronization, the case < q < 1 to partial synchro- 
nization, while for q = there is no synchronization at 
all in the system. 

In the N — * oo thermodynamic limit of the equilib- 
rium dynamics (t — * oo, so initial transient effects are 
lost) Kuramoto and Nishikava proved the existence of a 
K c critical coupling. For a Gaussian distribution of the 
oscillators natural frequencies, characterized by a D dis- 
persion they got m: 
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For K < K c the only possible solution gives q = (no 
synchronization) while for K > K c a stable solution with 
q ^ appears. Thus, the main result is that for a popula- 
tion of globally coupled non-identical oscillators a partial 
synchronization of the phases is possible whenever the 
interaction among oscillators exceeds a critical value. By 
adopting the Kuramoto model for biological populations 
the observed synchronization phenomenon are conceptu- 
ally understood. 



III. RHYTHMIC APPLAUSE VERSUS THE 
KURAMOTO MODEL: OPEN QUESTIONS 

As we argued in the introduction it is convenient to use 
the Kuramoto model in order to describe the dynamics 
of the rhythmic applause. If the Kuramoto model is ap- 
plicable to the phenomenon there could be two cases: i.) 
the average coupling between the spectators is smaller 
than the critical one and no synchronization appears, ii.) 
for larger coupling than the critical one synchronization 
would gradually evolve. The critical coupling (as dis- 
cussed in section 2) is governed by the dispersion of the 
spectators natural clapping frequencies (4). 

Studying carefully the rhythmic applause one might 
raise immediately many questions which are not obvi- 
ously answered within a simple application of the Ku- 
ramoto model: 

• Usually at the beginning of the applause there is a 
long "waiting" time without any synchronization, 
and with no increase in the order parameter. The 
partial synchronization is evolving after that sud- 
denly and achieves its maximal value in a short 
time. This should be not the case if we are in the 
K > K c limit. One would expect in this limit a 
continuous increase in the order parameter right 
from the beginning of the applause. 
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• Why synchronization already achieved is lost after 
a time, and why might reappear again? Loosing 
synchronization should not happen in the K > K c 
limit. 

In order to understand more deeply the phenomenon 
and to answer the questions which are not obvious within 
the framework of the Kuramoto model we first considered 
an experimental study. 

IV. EXPERIMENTAL STUDY 

Two main experiments were considered: 

I. The applause after many good theater and opera per- 
formances (in Romania and Hungary) was recorded, dig- 
itized and analyzed in several aspects. Recordings both 
by a microphone hanging form the ceiling of the concert 
hall and in the neighborhood of randomly selected indi- 
viduals were considered. 

II. Well-controlled clapping experiments were carried 
out on a group of 73 high-school students. We also inves- 
tigated the clapping frequencies of one individual during 
a one week period. 

A. Experimental method 

I. By digitizing the signal (Fig. la) we obtained a 
recording of fluctuating voltage with both positive and 
negative values. The zero level (as the mean-signal level) 
was determined and the square of the signal relative to 
this level was computed (Fig. lb). This signal would 
roughly correspond to the noise intensity variation, but 
due to the short sampling time it is definitely not the one 
that our human perception can follow. The average on a 
relatively short-time period (« 0.2s) was considered and 
we got the signal (Fig. lc) which describes well the useful 
noise-intensity variations during the applause. Already 
at a first-look one can easily detect the part where the 
rhythmic applause appears and the periodic noise inten- 
sity variation during synchronized clapping. 

This short-time averaged signal was analyzed in several 
aspects: 

1. A long-time averaged noise-intensity was computed 
by averaging the signal on a time interval of approx- 
imately 3s. 

2. An experimentally computable order-parameter, 
Qexp was calculated. This order-parameter was de- 
fined in a very analogous way with the q order pa- 
rameter (3) in the Kuramoto model. At each time- 
step q exp is calculated as the maximum of the nor- 
malized correlation between the s(t) signal and a 
harmonic function 



IXt s(t)sm(2n/T + 4>)dt\ 
q exp (t) = max {T , 0} <^ S 

(5) 

(in the above formula the values of <f> should span 
all possible initial phases between and 27r, and T 
should vary between two reasonably chosen limiting 
values, T min and T max .) Taking T min = 0.1s and 
T m ax — 5s q exp was numerically computed. 

3. Finally, for recordings taken in the neighborhood of 
one spectator one can observe that the short-time 
averaged noise intensity curves have some evident 
local maximums corresponding to the clapping of 
the individual. In this case we computed the time- 
interval between the clearly distinguishable groups 
of maxima as a function of time. This will charac- 
terize the clapping frequency of the chosen specta- 
tor. 

II. In our well-controlled clapping experiments we 
studied the distribution function for the natural clap- 
ping frequency of non- interacting students. Separated 
by the rest of the group, students were first encouraged 
to clap in the manner they would right after a good per- 
formance (mode I clapping). Second, they were asked to 
clap in the manner they would do it during synchronized 
clapping (mode II clapping). On one student we also per- 
formed investigations on the both two ways of clapping, 
sampling 100 times during a one-week period. 

B. Experimental Results 

I. A characteristic part of the global applause where 
the synchronization is formed and destroyed is presented 
as a function of time on Fig. 2a. On Fig. 2b. we present 
simultaneously the recording made in the vicinity of one 
spectator. The order parameter variation computed by 
equation || from the global signal is plotted on Fig. 2c. 
On Fig. 2d we show the variation of the long time moving 
average (window size: 3s) for the global noise intensity. 
Finally on Fig. 2e we present the computed clapping 
period as a function of time for the chosen individual. 
One can notice immediately that during the rhythmic 
applause (central part) the long-time averaged noise in- 
tensity has a clear minimum, the order parameter has a 
maximum (as expected) and the clapping period of an 
individual presents also a maximum. These results are 
stable and qualitatively consistent with all the 47 record- 
ings we have studied. 

II. For the normalized distribution of the natural clap- 
ping frequencies of the group of 73 students the results 
are given on Fig. 3. With continuous line we plotted 
the distribution obtained for mode I clapping and with 
dashed line the results for mode II clapping. Both distri- 
butions are roughly Gaussian. The maximum of the first 
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distribution is located at a higher frequency (approxima- 
tively two times higher) and the peak is wider than the 
one for mode II clapping. Calculating the dispersion for 
the two distributions, we got both the dispersion (D) and 
the relative dispersion (D r = D/u) for mode II. clapping 
smaller: 
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Computing the normalized distribution for the ratio of 
mode I to mode II clapping for each individual we got 
also a one-peek distribution centered around the value of 
2 (Fig. 4). 

The normalized distribution of mode I and mode II 
clapping obtained from 100 measurements on one indi- 
vidual (Fig. 5) shows a similar behavior with the one 
presented in Fig. 4. We plotted with continuous line the 
normalized distribution obtained for mode I clapping and 
with dashed line the one for mode II clapping. The dis- 
tribution of possible natural frequencies during mode II 
clapping is again sharper than the one for mode I clap- 
ping. Again, the location of the maxima for the distribu- 
tion of mode I clapping is roughly double than the one 
for mode II clapping. 

V. DISCUSSION 

From the controlled clapping experiments (Figs. 3-5) 
we conclude that spectators can clap with two very dis- 
tinct clapping modes. During rhythmic applause they 
shift from their original high frequency mode I clapping 
to the slow frequency, mode II clapping. The above pic- 
ture is totally supported by the clapping rhythm of one 
individual during the rhythmic applause (Fig. 2e). Dur- 
ing this period-doubling process the dispersion (width) of 
the natural clapping frequencies distribution is decreased 
roughly to half. This is a very important step in un- 
derstanding the phenomenon from the viewpoint of the 
Kuramoto model. From the Kuramoto model we learned 
that the critical coupling, K Cl for synchronization is di- 
rectly proportional with the dispersion of the frequency 
distribution of the oscillators (4). When the dispersion 
is decreased the value of the critical coupling decreases 
too. When the clapping starts the spectators clap very 
enthusiastically and with high frequencies (mode I clap- 
ping). The frequency distribution is wide and the value 
of the critical coupling is large. The coupling between 
the spectators is lower than the critical value and syn- 
chronization is not possible. The spectators inconstiently 
realize this (or already inconstiently know the "game" of 
the rhythmic applause) and by roughly doubling their 
clapping periods they shift to mode II clapping. By this 
the dispersion of the natural frequencies of the specta- 
tors reduces to half and so does also the value of the 



K c critical coupling. The coupling among spectators can 
become in this way larger than K c and partial synchro- 
nization (rhythmic applause) appears. 

It is clear now why synchronization appears, but the 
question remains: why after a time this synchroniza- 
tion is lost? The answer is immediate by examining the 
long-time averaged noise intensity curve (Fig. 2c). One 
can realize that after synchronization occurs the average 
noise decreases and attains its minimal value. Enthu- 
siastic spectators are not satisfied with this and try to 
increase the average noise intensity level. This is possi- 
ble mainly by speeding up the rhythm, because the inten- 
sity of one clap cannot be increased over a given level any 
more simply by hitting stronger. Speeding up the rhythm 
increases the natural frequency distribution's dispersion 
and consequently increases by this the value of K c . The 
coupling among spectators becomes again smaller than 
the critical one and synchronization is lost. 

The above scenario is fully supported by the analyzed 
signal in Fig. 2. We learn from here that during the 
rhythmic applause the individuals in the crowd double 
their clapping period (Fig. 2e), and we observe how by 
this the order-parameter characterizing synchronization 
is increased (Fig. 2d). When the order-parameter is max- 
imum the average noise intensity (Fig. 2c) is minimum 
and spectators will increase this by increasing the fre- 
quency of clapping shifting again to mode II clapping. 

The clapping after a good performance is frustrated in 
some sense. When maximal synchronization is achieved 
the average noise intensity is minimum due to the slow 
frequency of mode II clapping. By increasing the clap- 
ping frequency and shifting to mode I clapping the value 
of critical coupling is increased and synchronization is 
lost. The spectators cannot achieve both maximal noise 
intensity and good synchronization within the same clap- 
ping mode, and this makes the system frustrated. 

A very simple computer simulation exercise on the Ku- 
ramoto model can visualize all the above results (Fig. 6). 
Considering an ensemble of N = 70 globally coupled ro- 
tators with a Gaussian distribution of their natural fre- 
quencies (uj = 2ns' 1 , D = 2-K/6.9S- 1 and K = 0.8s" 1 ) 
we follow up their dynamics. We associate a very small 
time length r = 0.01s and Iq = uj/uj intensity noise for 
each event when a rotator passes through an integer mul- 
tiple of phase 2ir. The global noise is obtained by simply 
adding all noises generated by the individual rotators. 
In the given setup we have K < K Cl and thus no syn- 
chronization should appear. However, at t\ = 21s time 
moment we double the natural oscillation period of each 
oscillators, and as observable in Fig. 6a synchronization 
gradually evolves. This leads to an increase in the order- 
parameter (Fig. 6b, computed now as in experiments af- 
ter equation |5|) , and a decrease in the long-time averaged 
noise intensity (Fig. 6c). Beginning with t<i = 35s we 
begin to linearly increase the rotators natural frequen- 
cies back to their original values. As observable, this 
will decrease again the order parameter and increase the 
averaged noise intensity. This computer exercise approxi- 
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mates well the experimental results from Fig. 2, and gives 
confidence in our analytical statements. 

Some interesting comments can be made at this point. 

1 . In many societies the phenomenon of rhythmic ap- 
plause is unknown. It seems that the "game" of 
the rhythmic applause has to be learned by each 
community. Spectators once inconstiently recog- 
nizing that lowering the clapping frequency allows 
synchronization and getting familiar to the phe- 
nomenon they will keep alive this habit. 

2. Usually, in huge open air concerts no clear rhyth- 
mic applause forms even after outstanding per- 
formances. This is mainly due to the small and 
non-global coupling existing in the system (the 
applause of far-away spectators is totally unde- 
tectable). Even by reducing to half the natural 
clapping frequencies dispersion it is still not pos- 
sible to achieve a coupling larger than the critical 
one. 

3. In communist times it was a common habit to 
applaud by rhythmic applause the "great" leader 
speech. During this rhythmic applause the syn- 
chronization was almost never lost. This is a very 
nice evidence of the fact that spectators were not 
enthusiastic enough and were satisfied with the ob- 
tained global noise intensity level, having no desire 
to increase it. Frustration was not present in this 
system. 

VI. CONCLUSIONS 

By recognizing two characteristic clapping modes of 
the individuals during the applause and by applying the 
results of the Kuramoto model the dynamics of the rhyth- 
mic applause is qualitatively understood. The rhythmic 
applause is formed by lowering (roughly to half) the nat- 
ural clapping frequency of each individual, and thus by 
lowering the width of the clapping frequency distribu- 
tion of the audience. The rhythmic applause is lost and 
might reaper again due to two main desires of the spec- 
tators which can not be both fulfilled at the same time: 
optimal synchronization and maximal applause in- 
tensity. The system is frustrated in this sense. When 
maximal synchronization is achieved the average noise 
intensity is minimum due to the slow clapping frequency. 
Increasing the average noise intensity by increasing the 
clapping frequency synchronization gets lost. The amaz- 
ing interplay between unsynchronized and synchronized 
clapping is a consequence of this frustration. 
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FIG. 1. Digitized global signal, voltage as a function of time (a). The processed signals: (6) and (c). On figure 6 we have 
the square of the signal from figure a, relative to the calculated average voltage level. On figure c we present the short-time 
moving average (window size: 0.2s) of the signal from figure b. 
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FIG. 3. Normalized distribution functions for mode I (continuous line) and mode II (dashed line) clapping frequencies on 
a sampling of 73 high-school students. 
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FIG. 4. Normalized distribution functions for the ratio of the frequencies of mode I and mode II clapping (sampling on 73 
high-school students). 
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FIG. 5. Normalized distribution functions for the frequencies of mode I (continuous line) and mode II (dashed line) clapping 
of one chosen individual (sampling 100 times during a one week period). 



'g 10 



8 



0.1 
0.0 



1 1 


1 i 


i 


i i 






^^^^ b. : 




i i i i i 


i 1 i 1 

c. 



15 25 35 45 55 65 



t(s) 

FIG. 6. Computer simulation of the Kuramoto model for N = 70 rotators (K = 0.8s _1 , uj — 2irs~ 1 and D — 27r/6.9s~ 1 ). 
We double the rotators periods at t\ — 21s and linearly increase the frequency back to the original value from ti — 35s. The 
noise-pulse given by one oscillator has Io = uj/lj intensity and r = 0.01s duration. 
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